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Abstract 

The Backlund transformations for the relativistic lattices of the Toda type and their 
discrete analogues can be obtained as the composition of two duality transformations. 
The condition of invariance under this composition allows to distinguish effectively the 
integrable cases. Iterations of the Backlund transformations can be described in the 
terms of nonrelativistic lattices of the Toda type. Several multifield generalizations 
are presented. 

1 Introduction 

In this paper we consider four classes of difference and differential-difference equations, 
which contain, in particular, the Toda lattice introduced by Toda the relativistic Toda 
lattice introduced by Ruijsenaars ||2|, their generalizations introduced by Yamilov 



and Suris 0] and their discretizations introduced by Hirota and Suris ||9|, |l^, 11]. 
The systematic account is given of the method proposed in the papers |l^, |l^ for studying 
of these classes of equations, and some classification results are presented. The main idea 
of the method can be described in a few words. 

For the given lattice equation e[q\ =0 we define, in a certain way, the pair of transfor- 
mations : q ^ Q and : q ^ Q. In general, variables Q and Q satisfy the distinct 
equations. We say that the lattice e[q] = admits the duality transformations T± if both 
variables Q, Q satisfy the same equation E[Q] = which is called the dual equation. 
Obviously, the composition of the duality transformations can be used for reproducing of 
solutions. Denoting iterations of the transformations Ti = r_r^7+ and Tj^^ = Tj^TZ^ by 
superscript I one obtains the commutative diagram displayed on the Figure 1. 

In other words, transformations Ti and define the Backlund transformations for the 
given lattice and its dual. This explains the connection between duality transformations 
and integrability since existence of the Backlund transformation is an indispensable feature 
of any integrable system. Quite analogously, the Backlund transformations for the KdV 
and mKdV equations are obtained by the composition of two slightly different Miura 
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Figure 1. Backlund transformation as tlie composition of duality transformations. 



maps and the Schlesinger transformation for PII is constructed from two substitutions 
into PXXXIV. 

The requirement that the variables Q and Q satisfy the same equation is stringent 
enough and ahows to distinguish effectively the integrable relativistic lattices. However, 
it does not work for the subclass of nonrelativistic lattices, which are characterized by 
the property = T_. In contrast to the relativistic case, the duality transformation is 
now irrelevant to integrability and cannot be used for classification. Nevertheless, this 
subclass does not require special treatment since it arises as a by-product of already 
obtained results. Namely, it turns out that the iterations of the transformation Ti for the 
relativistic lattice admitting duality transformation are described by some lattice of the 
Toda type and the same is true for their discrete analogues. 

The above scheme is applied to the discrete relativistic lattices in the Section ||. In the 
Section 2A we introduce the notions of the duality transformation and the dual equation, 
and prove that the dual equation belongs to the same class. In the Section ^ we study 
the Backlund transformation T; in more details and prove that it is equivalent to some 
discrete lattice of the Toda type. The classification of integrable equations is performed 
in the Section 2/1, several multifield generalizations are presented in the Section ^. The 
Section |3| devoted to the lattices of the relativistic Toda type is, in fact, exact continuous 
double of the Section ^. Joint of both discrete and continuous theories is performed in the 
Section 3.4. 



2 The lattices of the discrete relativistic Toda type 

This section deals with the difference equations on the real or complex variable q defined 
on the two-dimensional integer lattice. As a rule, we use abridged notation q = qmn, 
Qij = Qm+i,n+j- The shift operators on the first and second subscripts will be denoted Tm 
and Tn respectively. Next, it is convenient to introduce notation for the differences (see 
Figure 2) 

x = q-q-i,o, y = q-qo-i, z = q - q^i-i. 

Obviously, one of these differences can be expressed through the other two, e.g. z = 
X + U-ifi = y + 2;o,-i) and the following identity is valid 

(r„ - l)yo,i = {Tn - l)xi,o. (2.1) 

The uppercase letters Q,X,Y,Z (reserved for the dual variable and its differences) are 
used in quite similar manner. 
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Figure 2. Pattern for the discrete relativistic lattice. 

2.1 Duality transformations 

The lattices of the discrete relativistic Toda type 

(T^ - + (r„ - \)g{y) + (r„r„ - = (2.2) 

are the Euler equations for the Lagrangians of the form 

/: = ^(a(x) + 6(y) + c(z)) (2.3) 

m.,n 

where f = a' , g = b', h = c' . In this section we assume that f'g'h' 7^ in order to eliminate 
the case of the nonrelativistic lattices, so that equation (|2.2| ) involves 7 nodes of the lattice 
as shown on the Figure 2. It is clear that the roles of the shifts Tm,Tn and T^^T~^ in 
(|2.2| ) are equal (actually, one can assume them as the generators of the regular hexagonal 
lattice). 

Equation ( |2.2D can be rewritten in two equivalent forms of the momentum conservation 
law: 

(T„ - l)(/(x) + hM) + (T„ - l){g{y) + h{z)) = 
^ (T„ - l)(/(x) + h{z)) + (T„ - l){g{y) + /i(zi,o)) = 0. 

This allows to introduce the pair of transformations acting on the differences: 

T+: (X,y_i,o) =T(x,yo,i), T_ : (Xo,-i, F) = T(xi,o, y) (2.4) 

where the mapping T : (x, y) {X, Y) is given by the formulae 



X = g{y) + h{x + y), Y = -f{x) - h{x + y). 



(2.5) 
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Obviously, the variables Q and Q corresponding to X, Y and X, Y are defined only 
up to the addition of an arbitrary constant. However, the equations for these variables 
contain only differences and can be derived by solving (p.4|), ( |2.5| ) with respect to x, y and 
using the identity ( |2.lD . Generally, these equations are distinct from each other. 

Definition 1. The lattice ^2.^ ) admits the duality transformations /\2.^ , l \2.3^ if the map- 
ping T is invertible and both variables Q, Q satisfy the same lattice which is called dual to 

The following Theorem characterizes equations admitting duality transformations in 
terms of the mapping T. It also demonstrates that the duality transformations do not lead 
out off the class ( |2.2D . The immediate corollary is that the equation which is dual to the 
dual equation coincides with the original one. 



Theorem 1. Equation (2.<:) admits duality transformations if and only if the inverse of 
(2.5) is of the form 

x = GiY)+H{X + Y), y = -F{X)- H{X + Y). (2.6) 

In this case the dual equation is of the form 

(r„ - l)F{X) + {Tn - l)G(y) + {T^Tn - l)H{Z) = 0. (2.7) 

Proof. Let be of the stated form, then one can easily check that the elimination of 
x,y by means of the identity ( p.l| ) brings to equation ( |2.7| ) for both transformations ( |2.4D . 
Conversely, assume that equations for Q and Q coincide. Let be of the form 

x = cD(X,y), y = ^{X,Y) 



then inverses of the transformations ( 2A ) are given by the formulae (the tilde in the second 
one is omitted) 

T+i; x = a>(X,y_i,o), yo,i = ^(^,^~i,o), 
TZ^: xi,o = ^(Xo,_i,y), y = ^(Xo,_i,y). 

The identity ( |2.lD yields the equations 

«>(Xi,i, yo,i) - $(^1,0, Y) - ^(Xi,o, Y) + ^(X, y_i,o) = 0, 

^>(x, yo,i) - HXo,-uY) - ^'(Xi,o, yi,i) + ^{x, yo,i) = o 

which must be equivalent to each other. In order to compare these equations, rewrite the 
last one in the form 

^>(x, yo,i) - $(x + y_i,o - y, y) - ^(Xi,o, + yo,i - Xi,o) + ^{x, yo,i) = o, 

so that both equations contain the variables Xi^i, Xi^, X, Yq^i, Y, Y^i^q. Now let us 
consider Xi^i as function on the rest variables from this set, then 

dYo,MXi,uYo,i) 



dYo^i 9xi,i$(Xi,i,yo,i) 

^ Oyp., [^{x, yo,i) - ^(Xi,o, + yo,i - Xi,o) + ^{x, yo,i)] 
axii^(Xi,o,Xi,i + yo,i-Xi,o) 
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The last equality must be satisfied identically. Differentiating it with respect to X and 
Xi^o yields 

5x^o,i(^(^.^0,l) + ^(^,^^0,l)) =0, 9xi,o5yo,i^(^l,0,^l,l+>0,l-^l,o) =0 
and hence 

^{X,Y) + '^{X,Y) = G{Y)-F{X), '^{X,Y) = K{X) - H{X + Y). (2.8) 
Analogously, differentiating the relation 

ay-1,0 _ dx^{x,Y^i,o) 

_ dxMx, yo,i) - ^x + y_i,o -y,y) + ^{x, Yq,,)] 

dy_,^,^X + Y_i,o-Y,Y) 

with respect to Y yields 

dxdY^X + Y_i^o-y,y) = <^iX,Y) = L{Y) + M{X + Y). 
Comparing with ( |2.§| ) completes the proof. I 



2.2 Classification theorem 

The Definition |^ turns out to be severe enough and allows to obtain the finite list of 
integrable equations ( |2.2D . In virtue of the Theorem |^ it is sufficient to find all functions 
f,g,h such that inverse of the transformation ( p.5|) is given by (|2.6| ). This means that 
the Jacobian A = f'g' + g'h' + h' f of the map (|2.5D must be nonzero and the following 
identities must hold 

XXY + yXY = 0, XXY = XXX, UXY = UYY ■ 
These three relations are equivalent. Indeed, the Jacobi matrix is 

XX yx\^l( -h' f + h' 

XY VY ) A V -g' - h' h' 

that is XX = —Uy- Straightforward computation proves that functions /(x), fl'(y), h{x + y) 
must satisfy the equation 

f" a" h" 

i9' + h')^ + {f' + h'f- = {f' + g')^. 

The designations f = 1/u, g' = 1/v, h' = 1/w rewrite it in more convenient form 

[v{y) + w{x + y)]u'{x) + [u{x) + w{x + y)]v'{y) = [u{x) + v{y)]w'{x + y). (2.9) 

The classification problem is reduced to solving of this functional equation. Of course, 
functions n, v, w can be multiplied by an arbitrary constant and the linear transformation 

x = c{x - xo), y = c{y- yo) <^ Qm^n = c{qm,n " mxQ - nyo - const) 
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can be applied, as well as the permutation of the x, y, z axes. 

At first let us consider some degenerate cases. Assume that two of three functions are 
constant, say u and v. Then ( |2.9| ) yields that either w is constant as well, or n = —v and 
w is arbitrary. In the first case the equation ( |2.2D is linear and in the second one it is of 
the form 

aiqifi + q-ifi - qo,i - qo~i) + Hqi,i -q) - Hq - q-i-i) = o 

and admits "integration": 

(^{qm+l,n — qm,n+l) + ^(^m+l.n+l — qm,n) = Cra-n- 

Other degenerate case corresponds to the vanishing of the Jacobian, what is equivalent to 
u + V + w = Q and implies that all three functions are linear. In this case one can prove 
that equation ( |2.2D can be reduced to the equation on the variable p = yo,i/x. Further on 
we will not consider these degenerate cases. 

Lemma 2. The functions u, v, w satisfy the equations 

{u'f = Su"^ + 2au + e, iv'f = 5v'^ + 2Pv + e, iw'f = 5w'^ + 2'yw + £. (2.10) 
Proof. At first prove that functions u{x) and v{y) satisfy the equation 

(u" -v"){u + v)-{u')'^ + {v'f = k{u-v), A: = const. (2.11) 



Let us eliminate w from ( [2.91 ). Applying the operator dx — dy one obtains the linear system 



on w, w' : 



u' + v' —u — v\fw + u\ , s f u' 
n II I I W I \=[u-v) 
U — V V — u J \ w J 

Its determinant A is exactly the left hand side of the equation ( |2.1l| ). If it is identically 
zero then ( |2.11| ) is proved, otherwise one finds 

w + u= ^^^^^iuu" - (u'f + u"v + u'v'), w' = ^^-—^(u"v' + u'v") 
A A 

and consequently {{u — v)/ A)y{uu" — {u')'^ + u"v + u'v') = 0. Assume that the expression 
in the second bracket vanishes. If u' ^ then v' = —u"v/u' + u' — uu" /u' , v" = —u"v'/u', 
but then, as one easily checks, A = 0. If = then w + u = 0, that is we come to the 
degenerate solution excluded above. Therefore ((u — v)/A)y = 0. Due to the symmetry 
between u and v one can prove analogously ((n — v)/ A)x = and obtain ( 2.11| ). 
Further on, rewriting ( |2.1l| ) in the form 

u' \ _v" + k {v'f + 2kv 
^u + vj^ u + v {u + vy ' 

multiplying by u' / (u + v) and integrating with respect to x yield 

{u'f = 6iy)iu + vf - 2{v" + k){u + v) + {v'f + 2kv. 
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Replacing u and v one obtains 

[v'f = ~6{x){u + vf - 2{u" + k){u + v) + {u'f + 2ku. 



Subtracting one equation from another and using ( |2.11 ) one obtains 6 = 5 = const . 



Summing and dividing hy u + v give u" + v" = 5{u + v) — k. The separation of the variables 
yields (u')^ = Sv?' + 2au + £, {v')"^ = dv"^ + 2l3v + e, where a + /? = —k, and substitution into 
(2.11) proves e = e. The last of the equations ( 2.1C| ) is obtained in virtue of the symmetry 



of X, y, z axes. ■ 

It is clear that the solutions of equations ( p.lOD must satisfy also some additional 
relations. However their analysis is not in principle difficult, and the direct examination 
of all solutions brings to the following list. 



Theorem 3. The equations (2.i) admitting duality transformations are exhausted, up to 
the changes qm,n = c{qm,n — mxQ — nyo) and permutations of x, y, z axes, by the following 
sets of the functions f,g,h. In formulae (A), (B), (C) the parameters are constrained hy 
relation X + ^ + v = Q, and in (I) hy relation Xfiu = — 1. 



(A) 


f = 


X ' 




g 


V 


h = 


X 

z ' 


(B) 


f = 


/i coth X, 




9 


= i^cothy, 


h = 


A coth z, 


(C) 


f = 


2 ' 




9 


2 y-u' 


h = 


1 log 


(D) 


f = 


logx, 




9 


= logy, 


h = 


log(l - 1/z) 


(E) 


f = 


-e^ - 1, 




9 


= e-^ 


h = 


1 


(F) 


f = 


log(e^' — 1) 




9 


= \og{ey - 1), 


h = 


— log(e^ — 1 


(G) 


f = 


— log(e~^ - 


-1), 


9 


= \og[ey - 1), 


h = 




(H) 


f = 


log(A-He^' 


+ 1)), 


9 


= log(e-S'-l), 


h = 


loe ''"+^ 

log ^^'+^ 


(I) 


f = 


loe 




9 


= log ^""+1 


h = 



The duality transformations ([2.^), I{2.9^ link together the equations corresponding to so- 
lutions (B) and (C), (D) and (E), (F) and (G), while the equations corresponding to 
solutions (A), (H) and (I) are self-dual. I 

Notice, that the cases (A) and (B) are connected by the point transformation q = 
exp(2g). It is explained by the fact that the Lagrangian ^(//logx + i^logy + Alogz) of 
the equation (|2.2| ), (A) is invariant under the dilations q ^ Cq as well as under the shifts 
q ^ q + C. On the other hand, the inversions g — > g/(l — Cq) preserve the Lagrangian as 
well but the change q = 1/q which maps this group into the shift group does not bring to 
a new equation. 

2.3 The lattices of the discrete Toda type 



The method proposed in the Section |2.1| does not work for the lattices of the discrete 
Toda type which correspond to the case h = (this is equivalent to f'g'h' = 0, without 
loss of generality). Indeed, in this case transformations r+ and r_ coincide for arbitrary 
/, g and classification becomes impossible. However, we can dispense with it since these 
lattices arise as a by-product of already obtained results for the discrete relativistic lat- 



tices. Namely, we will demonstrate, using only few basic formulae from the Section 2.1 
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and without any comphcated calculations, that the iterations of the Backlund transfor- 
mation Ti = TZ^Tj^ are described by some discrete lattice of the Toda type. Hence we 
automatically obtain some list of integrable lattices, see Theorem |^ below. Probably, this 
list is exhaustive (cf. IC, 11, |l^), but unfortunately I do not know any classification 

|, which can be compared with this list. 



results, like Yamilov's Theorem 



Let us consider some lattice (2^) admitting duality transformations and denote itera- 



tions of Ti by superscript /, in such a way that tilde in the formula ( |2.4D corresponds to the 



value / — 1. It is possible to rewrite equations (|2.2| ), (2.7) in terms of the mixed variables 



Indeed, one obtains directly from ( |2.4| ), (^) the relations 

-'^o -1 = g{y) + h{xo-i + y), X''^ = g{yo,i) + h{xi^i + yo,i) 

and therefore equation ( |2.2p is equivalent to 

(r„-i)/(x) + x-i-Xo,_i = o. 

Analogously, in virtue of (|2.6|) equation ( |2.7D is equivalent to 
(r^-l)F(X)+xi,i-xl,o = 0. 

"Integrating" this equation with respect to m (recall that xi^o = (^m — 1)?) one obtains 
X = ip{{Ti — Tn)q + c) where function if is inverse of F. The constant c does not depend on 
m, but may depend on l,n, and it can be set to zero without lost of generality by means 
of appropriate shift of the variables q. Then eliminating X from the previous equation 
brings to the lattice of the discrete Toda type 

(T^ - l)/(x) - {TiT~' - - T~^Tn)q) = 0. (2.12) 

So, we have already proved that the Backlund transformation for the discrete relativistic 
lattice is governed by some discrete nonrelativistic lattice. However, the complete picture 
is even more rich: it turns out that the 3-dimensional lattice generated by the shifts T/, 
Tm, Tn contains 3 instances of the discrete nonrelativistic lattices and 4 instances of the 
discrete relativistic lattices. In order to see this let us remind that the roles of all shifts in 



equation (2^) are equal and therefore we can repeat the above calculation starting from 
the other set of mixed variables. More precisely, rewriting equations ( |2.2| ) and ( |2.7| ) in 
terms of y, Y or z, Z one obtains the equations 

(r„ - l)g{y) + Yzlo - y = 0, (r„ - l)G(y) + yo,i - yli = 0, 
(r„r„ - l)h{z) +Z-Z-^ = 0, {TmTn - 1)H{Z) + zl, - = 

and further elimination of Y and Z yields equations 

{Tn - l)g{y) - {TiTm - 1)V'((1 - T^^TZ,^)q) = 0, (2.13) 
{T^Tn - l)h{z) + {Ti - l)r?((r-i - l)q) = (2.14) 

where functions ^ and r/ are inverses of G and H respectively. 

The placement of the variables involved in equations ( 2.12] ), ( 2.13| ) and ( |2.14| ) is shown 



at the Figure 3. At this picture the variables q involved in the original equation (2.2) 
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n 




m 



Figure 3. Each of three discrete nonrelativistic lattices involves 2 pairs of the same shape (black 
and white). Each of four discrete relativistic lattices involves 3 distinct pairs, number of white 
pairs is even (0 or 2). 



lie in the plane (mn) and are marked by black symbols, while the variables obtained by 
means of the Backlund transformations are marked by the white ones. Equation ( ^l2|) 
constraints two black and two white circles (and, of course, the centre of the cube), squares 
and triangles correspond to ( ^.13 ) and ( 2.141) . 

Next, subtracting from the equation ( |2.2| ) two of three equations ( 2.12| ), (2.13) and 



(|2.14| ) we obtain some discrete relativistic lattice again, for example subtracting of ( p. 131) , 



( 2.14 ) yields equation 



(T^ - l)/(x) - {Ti - l)r?((r-i - l)q) + (T^T^ - 1)V^((1 







which involves variables in the plane {Im). Two other choices bring to the discrete rela- 
tivistic lattices in the plane {In) and in the plane spanned over the vectors (0, 1, 1) and 
(1,1,0). 

In conclusion of this section we present the list of the discrete lattices of the Toda type 
which is obtained by direct examination of the equations listed in the Theorem 0. 



Theorem 4. The equations { 2.1i ), ( 2.1^ , ^-14 ) o,'^^ equivalent, up to the renaming of 
the shifts and the linear changes qm,n = c{qm,n — mxQ — nyo), to the lattices 



1)- = {Tn 

X 



y 
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(T^ - l)e^- = (r„ - l)e^ 

(r^-i)iogx = (r„-i)iogy, 

(r„, - 1) log(l - l/x) = {Tn - 1) l0g(l - 1/y), 

(r„, - 1) log(e- - 1) = (T„ - 1) log(e2' - 1), 
{T^-l)x = {Tn-l)\oz{e^ -I), 

^^--l)^°s(J^)=(T„-l)log(^ 

3 The lattices of the relativistic Toda type 

Now let us consider the differential-difference equations on the variable qn{t). As before 
we will omit the subscript n : q = Qn, qi = qn+i and use uppercase letters for the dual 
variables. Instead of differences x,y,z we consider the quantities 

p = Q, y = q-q-i 

which satisfy identity 

DtVi = {Tn - l)p. (3.1) 

3.1 Duality transformations 

The lattices of the relativistic Toda type 

p = r{p){h{yi)pi - h{y)p--i + g{yi) - g{y)), y = p-p_i (3.2) 
are the Euler equations for the Lagrangians of the form 

C= f dtY,{a{p) - h{y) - c{y)p), (3.3) 

where r = l/a" , g = b' , h = c' . In this section we assume that the nondegeneracy condition 
/i / 0, l^'l + / is fulfilled. 

The term c{y)p of the Lagrangian is equivalent, up to the total divergence, to c{yi)p 
and this results in two equivalent forms of the momentum conservation law: 

Dtia'ip) - c(yi)) = (r„ - l)(/i(y)p-i + giy)) 
^ Dt{a'{p) - c{y)) = (r„ - \){h{y)p + g{y)). 

This allows to define the pair of transformations 

T+: (P,y)=T(p,2/i), r_: {P.^.Y) = T(jp,y) (3.4) 

where mapping T : {p, y) — > (P, Y) is defined by the formula 

P = h{y)p + g{y), Y = a' {p) - c{y) . (3.5) 



44 



V.E. Adler 



Notice, that the lattices ( |3.2D are the continuous hmit of the discrete relativistic lattices 
(|2.2| ). Indeed, let us consider the family of the functionals of the form (|2.3| ) 

l^e = ^(ea(xm,„/e) - eb{ym,n) + k{ym,n) - k{Zm,n)), 

m.n 

and let k' = c, qm,n = Qn{t), t = me. It is easy to see that passage to the limit e ^ 
brings exactly to the Lagrangian (^]^). Transformations (^^) also are obtained from ( |2.4D 
by this limit. 

Definition 2. The lattice l \3.^ ) admits the duality transformations l \3.^ , l \3.3^ if the map- 
ping T is invertible and both variables Q, Q satisfy the same lattice which is called dual to 



Notice that functions a, 6, c are defined up to some linear transformations and therefore 
mapping ( |3.5| ) is not uniquely defined. However, it is easy to see that this arbitrariness 
corresponds to the linear changes of the variables P, Y and all lattices dual to the given 
equation ( |3.2D are equivalent under the changes Qn aQn + f3n -\- jt. 



Theorem 5. The lattice {3.i) admits duality transformations if and only if the inverse 



of the mapping ([3.^ is of the form 

p = H{Y)P + G{Y), y = A'{P)-C{Y), H = C . (3.6) 
In this case the dual equation is of the form (R = 'i/A") 

P = R{P){H{Yi)Pi-H{Y)P.i + G{Yi)-G{Y)), Y = P - P_i. (3.7) 
Proof. Let T~^ be of the form (|3.6| ). Then elimination of p, y from ( ^.4| ) by means of 



identity ( pj| ) brings to equation (3/7) in both cases. 

In order to prove inverse statement we need to compare equations on the variables Q 
and Q. Straightforward calculation proves that the variables P,Y satisfy the equations of 
the form 

P = r{p){h{yi){P-P^i)+A{p,yi){pi-p)), Y = P - P_i, (3.8) 

where A(p, y) = h'^{y) + :^{h'{y)p + g'{y)) and yi have to be expressed in terms of 
Pi,P,Yi,Y. Therefore, the dual equation must be linear in P-i. Analogously, considering 
equation on Q one proves that the dual equation must be linear in Pi. Since in (|3.8| ) only 
pi depends on Pi, hence mapping T must satisfy the condition d'^p/dP^ = 0. This proves 
the first formula in (^.6|). The second one follows from the property dy/dY = —dp/dP 
which is evident from the structure of the Jacobi matrix 

/ dp/dP dy/dP \ ^ 1 / h{y) l/r{p) \ 

\dp/dY dy/dY ) A{p,y)\h'{y)p + g'{y) -h{y)J ^''•^^ 

where A 7^ by Definition |2[ ■ 

As in the discrete case we see that the original lattice is dual for its dual and the 
composition TZ^T^ defines the Backlund transformation for (|3.2|) (see Figure 2). We will 



study this Backlund transformation in the Section 3.3 
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3.2 Classification theorem 

Theorem 6. The lattices l \3. ^ ) admitting the duality transformations are characterized 
by the following equations for the coefficients: 

r = r2p'^ + rip + rQ, 

2 2 / 2 (3-10) 

g = r2g + Rig + Rq — roh , h = 2r2gh — rih + Rih. 



The coefficients of the dual lattice satisfy equations 

(3.11) 



R = raP^ + RiP + Rq, 



G' = raG^ + riG + tq - RqH^ , H' = 2r2GH - RiH^ + riH. 

Proof. Straightforward calculation proves that necessary and sufficient condition 
d'^p/dP^ = is equivalent to the following relation involving the functions r{p),g{y), h{y) : 

h(h"p + g") - h'Qi'p + g') + 2rh^h' - r'h^{h'p + g') = 0. (3.12) 

The second derivative of ( p. 12 ) with respect to p is r"'{h'p + g') = 0. By the condition of 
nondegeneracy, this implies that r is a polynomial and its degree is less than 3. Dividing 
( 3.12| ) by /i^ and integrating the result with respect to y one obtains 

h'p + g + 2rh^ — r'h{hp + g) = h{a + Rip) 

where a and Ri are some constants. Let r = rap^ + rip+rg, then collecting the coefhcients 
of p in this relation gives the system 

g' = ah + righ — 2roh'^, h' = 2r2gh — rih'^ + Rih 

which is equivalent to ( |3.10 ) modulo common first integral 

^25^ + Rig + Ro — figh + tq/i^ = ah. (3.13) 

The formula for R can be easily obtained from the relation R = A{p,y)r. Thus, equa- 
tions for g and h are uniquely defined by the coefficients of the polynomials r and R. Since 
due to the Theorem ^ the duality relation is symmetric, hence equations for G and H can 
be written automatically. I 

Remark. The question about value of the first integral for the system ( |3.1l| ) is a bit 
more complicated. However, this value can be calculated by sequential elimination of g, h 
and p from the formula ( 3.13| ) using relations g = P — hp, h = RH/r and p = HP + G. 
This brings to relation 

r2G^ + riG + ro - RiGH + RqH^ = aH, 

that is, the value of the integration constant a for the dual lattice is the same. Notice, 
that actually role of this constant is very important since it may depend on n. This brings 
to the integrable lattices containing an arbitrary parameter in each node ||l^ (it plays role 
of the discrete spectrum when constructing solitons). For sake of simplicity we will not 
consider this generalization. 
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By use of the first integral ( |3.13|) the system (l3.1C| ) is reduced to the equation 

{h'f = {rl - 4r2ro)/i^ + (4ar2 - 2Riri)h^ + (iif _ 4r2Ro)h'^ (3.14) 

which is solved in elementary functions. Consideration of all the possible choices of pa- 
rameters and the branches of solutions brings to the following result. 

Theorem 7. Up to the transformations aqn + f3t + 'jn, t — > 6t, the nondegenerate 

lattices admitting the duality transformations are exhausted by the list {y = p—p-i) : 



(a) 


p 


= piey^ — p-iey 


{b) 


p 


( Pi p-1 
= p[ 

\yi y 




p 


_ f Pi 


id) 


p 


= p{p + l)l — 



+ yi-y ^ 

+ u(e 



p = p{p- fi) 



1 + fie-y 

P-1 
y 

pi p~i 



fi + ey^ ^i + ey 



if\ ■ / 2 , . Pi-yi p-1 -2/ 

{f^.) P={P +P)\ ,,^^,2 - „ , o 



{9^l) P=^(p^ + 1-^ 



1^ 2 , 1 2\/^Pi~s™hyi sinhy 



^ + cosh yi fi + cosh y 



The duality transformations ^.4 ); (3-5) link together equations (a) and (b), (d) and (eo), 
(f/j.) for /U 7^ and (g±i), while the rest equations are self-dual. ■ 

3.3 The lattices of the Toda type 

As in the discrete case, the subclass of the Toda type lattices 

p = r{p){f{yi) - f{y)), y = p - p^i (3.15) 

must be considered separately, since transformations and T_ coincide. It should be 
mentioned that in this case the duality transformation was introduced by Toda [||. It is 
given by the formula 

P = f{yi), Y = a'{p), a" = l/r 

and the coefficients of the dual lattice are defined by the formulae /' = R{f), F{a'{p)) = 
p. For example, the Toda lattice q = exp{qi — q) — exp((7 — q-i) is dual to the lattice 
Q = Q{Qi-2Q -Q^i). 

In contrast to the relativistic case, the duality transformation is irrelevant to integra- 



bility and cannot be used for the classification of the lattices ( 3.15 ). For the first time this 



problem was solved by Yamilov in the framework of the symmetry approach. 



Theorem 8 (Yamilov, [Q]). The Toda type lattice l ^-l^ admits the higher symmetries 
^ff 

r{p)=r2P^ + rip + ro, f = r2f^ + R^f + R^. (3.16) 
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This result demonstrates that integrable lattices ( ^.15 ) can be obtained from the in- 



tegrable lattices of the relativistic Toda type by passing to the limit /i = in equations 
( 3.10| ). More interesting link between this two classes of equations can be established along 



the same arguments as in the Section 2.3 



Let us denote iterations of transformation Ti = TZ^Tj^ by the superscript / and let tilde 
in the formula (|3.4D corresponds to the value I — 1. Then formulae (|3.4D , (|3.5|) , (|3.6| ) take 
form 

P_i = hiy)p^i + giy), = h{y{)pi + g{yi) (3.17) 

p^ = H{Y)P_i + G{Y), p = H{Y)P + G{Y) (3.18) 



and therefore equations ( |3.2| ), ( |3.7D can be rewritten as a coupled lattice of the Volterra 
type 

p = r(p)(p-i -P_i), P = R{P){pi-p'). 
Assuming 

P = f{qi-q^). f = R{f) (3.19) 
we obtain the Toda type lattice 

q = r{q){f{qi' - q) - fiq - q-i)) (3.20) 

for the variables situated along the line I + n = const . Since functions r, R were described 
in the Section ^, we immediately repeat the Yamilov's result ( p.l6| ). 

3.4 Nonlinear superposition principle 

Recall, that in the discrete case the 3-dimensional lattice generated by the shifts T^, T^, 
Tn contains 4 instances of the discrete relativistic lattices and 3 instances of nonrelativistic 
ones (see Figure 3). In the continuous case the picture is more bare: the variables qin form 
the 2-dimensional lattice containing 2 instances of the relativistics lattices, 1 instance of 
the Toda type lattices and 1 instance of the discrete Toda type lattices, as shown on the 
Figure 4. 

The shift T„ corresponds to the lattice (3^) and the diagonal shift T^^^T„ corresponds 



to the lattice ( 3.20D , as was proved in the previous section. The following Theorem demon- 



strates that the condition of the commutativity of these shifts is equivalent to some discrete 
lattice of the Toda type which can be interpreted as the nonlinear superposition principle 
of the equations ( |3.2| ) and ( p. 20 ). 



Theorem 9. The variables q,q^^,q±i are related by equation of the form 

{Ti - l)c(g-i -q + 5) + {Tn- l)c(g - q^i) = 0. (3.21) 

Proof. The function inverse to A' satisfies equation ip' = R{f), that is f{y) = ip{y + e). 
Therefore 

CiY) = A'{P) - 2/1 = A'(/(gi - q^)) -q,+q = q-q^+s, 
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Figure 4. Equations associated with the relativistic lattice. 

that \s Y = s{q — + e) where s is the inverse function to C. It is easy to prove that s' 
satisfies the same equation ( |3.14| ) as h = c' , that is s{y) = c{y + e) + const . 
Next, ehminating p from the formulae 

Y = a'{p)-c{yi), Y-^ = a'{p)-c{y) 

brings to Y~^ — Y = (T„ — l)c(y) and this completes the proof. ■ 

It is clear that parameter 5 can be set to zero by the proper choice of the function / in 
formula ( ^3.19 ) or, equivalently, by the shift — > q\^ + 51. Therefore, the terms in ( ^.21 ) 
are symmetric and the shifts T; and r„ plays the equal roles, although their origin was 
different. This suggests that T/ corresponds to some relativistic lattice as well. Indeed, 
eliminating P from the equation p = r[p){P~^ — P-i) by means of the relations ( p. 18 ) 
brings to equation 



p = r{p) 



p. 



+ 



H{Y) H{Y 



+ 



G{Y) G{Y- 



H{Y) H{Y-^) 



and since Y = c{q — + 5) we obtain a relativistic lattice again. The coefficients of this 
equation coincide with /i, g up to the shift of argument. 



3.5 Example: Heisenberg chain 

In the previous sections the discrete symmetries of integrable lattices were studied. Now 
we will briefly discuss the continuous symmetries. I restrict myself by example of the 
equation (fo)[|- This choice is motivated by the link between (fg) and the Heisenberg 
chain which became very popular due to its applications in discrete geometry 17, 18]. 
Remind that this model reads 



St 



as X 



Si 



+ 



S-l 



l + (s,si) l + (s,s_i) 



+ h 



Sl+ s 



S + S-l 



l + (s,si) l + (s,s_i) 



(3.22) 



^These results were obtained in collaboration with Professors A.B. Shabat and A. P. Veselov. 
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where s G M'^, {s, s) = !,(,) and x denote standard scalar and vector products respectively, 
a and b are arbitrary constants. Up to the author knowledge this equation was introduced 
by Sklyanin [^] in the case 5 = and by Ragnisco and Santini |2C] in the general case. 
The continuous limit of the Sklyanin lattice is the Heisenberg model 

st = sxsxx, {s,s) = l. 

It was noticed that r-matrices for these models coincide and, more generally, this property 
can be accepted as a definition of correct discretization for a given continuous equation 

H- 

In order to obtain (fo) let us consider the complexification s € M'^ ^ s G and the 
stereographic projection 

s = S{u, v) = (1 — uv, i + iuv, u + v). (3.23) 



It is convenient to represent the flow ( |3.22| ) for arbitrary set of parameters a, 6 as a linear 
combination of the flows corresponding to the sets a = z, 5 = ±1. These flows are given by 
the following formulae in terms of the variables u, v : 

(^1 -u){u-v) {u - v){v - V-l) 
ut+ = , vt^ = , (3.24) 

Ul — V u — «_l 

(n-1 -u){u-v) {u-v){v - vi) 
ut^ = , vt_ = . (3.25) 

U-l — V u — Vl 

The lattice ( p.24| ) appeared in |22, ^] for the first time. 

Next, notice that elimination of P in virtue of ( |3.19D brings the formulae ( 3.17 ) to the 
form 

. ^ /(gi -g^) -g) .1 ^ fiQ-Q-i) -ajq^ -Q-i) 

h{qi-q) ' /i(g^-gii) 

It is easy to see that this system corresponding to the case (fo) (/(y) = g{y) = — 
h{y) = 1/y^) coincide with ( ^.24 ) if we assume q = u, = v and t = t+. Since the lattice 
( 3.25| ) is obtained by reflection n —n we immediately come to the following statement. 

Proposition 10. Both variables u and v satisfy the lattices 



,2 



qt±t± - qt± 7- -y2 - Tz — : — ^2 ~ z 1 + 



^ V (g±i - qY {q - q^fiY q±i-q ' g - g^i . 

in virtue of the equations i\3.24) , fi3.2$ ). 

The lattices (II2I , ( ^) are Hamiltonian. Their Hamiltonians are 

H+=y log , i7_ = > log 

^-^ u — V ^-^ u — V 

n n 

respectively and the Poisson brackets are of the form 

{■Um, Vn] = {Un - VnfSmn-, {Um,Un} = {Vm, Vn} = 0. (3.26) 

The next proposition can be proved by straightforward calculations. It demonstrates 
that the lattices ( 3.24| ), ( p.25| ) are the symmetries of each other and the shift (n„, v„) 
{un+i,Vn+i) defines the Backlund transformation for some hyperbolic system. 
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Proposition 11. The Hamiltonians Hj^ and i7_ are in involution, the vector fields dt^ 
and dt_ commute and the variables Un,Vn satisfy the system 

^ut+ut_ 2vt^vt_ 
ut+t_ = ut+ - ut_, vt^t^ = h vt+ + vt_. (3.27) 

U — V V — u 

In terms of the vector ( 3. 23^ ) this system reads 

st+t- + {st+,st_)s + is X {st+ + st_) = 0, {s,s) = 1. 

Commutativity of the flows dt_^_, dt_ allows to construct compatible zero curvature 
representations 

Wt+ = U^W - WU+, Wt_ = U^W - WU-, 

so that tr Wn ■ ■ ■ Wi generates the common first integrals of the both lattices under the 
periodic boundary conditions un = uq, vn = vq. The matrices W,U^ are given by the 
formulae 

W = \I + P{u, v), U+ = ^{P{u, v-i) - P{v-i, u)), 



u- 



1 



2(A + 1 

where P denotes the projector 



iPiv,u^i)-Piu^i,v)) 



P{u,v) 



1 



-V uv 



li — u \ — 1 u 

The Poisson structure ( |3.26| ) can be written in the r-matrix form 
with the same r-matrix as for the Heisenberg model: 



/ 1 \ 

1 10 

A-/i 10 

\ 1 / 



Equation ( |3.27| ) gives an example of hyperbolic symmetry. Any lattice (^]^) of the rela- 
tivistic Toda type admitting duality transformations possesses also evolution symmetries, 
the simplest representative is given in the following Theorem. 



Theorem 12 ( |12[| ). The Lagrangian (3^) admits the variational symmetry of the form 
Qr = r{p){h{yi)pi + h{y)p^i + g{yi) + g{y)) + Rip^ (3.28) 



if and only if the coefficients r,g,h satisfy the system ( 3.1L ) 
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Obviously, one can rewrite this symmetry as the partial differential equation if use the 
lattice ( ^^ ) itself and assume q = u, qi = v : 

Ut + utt = 2r{ut){h{v - u)vt + g{v - u)) + Riuj, 

Vr - vtt = 2r{vt){h{v - u)ut + g{v - u)) + Rivl- 

Returning to the system (|3.24| ) one obtains the symmetry (t = 



Utt = 



vt 



1 



iVr + Vtt = 2Vt 



{u — vY 

Ut 

{u — vy 



u — V 

1 

U — V 



(3.29) 



which is equivalent, in the geometrical terms ( |3.23| ), to the modified Heisenberg model ||2 



S X Stt - ■^st{st,st), 



{s,s 



1. 



(3.30) 



So we have proved that the flows defined by the lattice ( p.22 ) at a = i,b = 1 and equation 
( 3.30| ) commute. One can reformulate this result as follows: 



Proposition 13. The shift s — > si in the Heisenberg chain {3.2i) for a = i,b = 1 defines 
the t-part of the Backlund transformation for equation ( \3. 3(\) . 



Notice, that according to (3.24) this shift is equivalent to solving of Riccati equation 
on the variable vi : 

{ui-u){u-v) {ui - vi){vi - v) 

Ut = , wi,t = - 



Ul — V ' Ul 

Another choice of dependent variables allows to rewrite ( |3.29| ) in the form 



2u^ 



lUl- 



Ul,tt 



l,t 



-Vtt 



2.2 



Ul — V Ul — V 

This system is equivalent to the Heisenberg model ar = a x att, (o", cr) = 1 in terms of the 



vector a = S{ui,v), which is related with the vectors (3.23) by the formula 
si + s — isi X s 



a 



l + {si,s) 



4 Multifield examples 

In the previous sections we assumed that the variable q was scalar {q G F, F = M, C). The 
brief analysis proves that the Theorems ^ and |^ remain valid for the vector case q S F^ 
as well. Thus, the multifield generalizations can be obtained by finding of the mappings 
T : F^^ with the special structure described in these Theorems. This problem is 

much more difficult than in the scalar case and the classification of such mappings is far 
from completeness. 

However it is not difficult to find some particular examples. The most simple are the 
generalizations of the discrete Heisenberg equation (A) and the Heisenberg lattice (fo) 
related to the Jordan triple systems. It was Svinolupov who recognized the role of the 
Jordan algebraic structures in the theory of integrable systems for the first time. I provide 
only necessary information about the Jordan triple systems, the interested reader can find 



more in |25, 27, 28, 29] and references therein. 
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4.1 Jordan triple systems 

The ternary algebra J with multiphcation {} : J'^ ^ J is cahed Jordan triple system if 
the following identities hold 



Some consequences are: 

{ab{aca}} = {a{hac\a} = {{aba}ca}, {a{6a6}c} = {{a6o}6c}, 

2{{ahc}hd} = {a{hcb}d} + {c{hab}d}. ^ 

Calculations are simplified by use of the linear operators Lab, ^ah-, ■ J ^ J defined for 
arbitrary elements a,b J as follows: 

Labic) = {abc}, Mab{c) = {acb}, Ma{c) = Maa{c) = {aco}. 

Notice, that identity ( |4.2| ) is equivalent to 

[Lab, Led] = L{cba}d - Lc{bad}- (4-4) 

We are especially interested in rational expressions. They are build from the inverse 
elements which are defined as = Ma^{a). In some Jordan triple systems the operator 
Ma is degenerate. In such cases the notion of inverse element can be partly substituted 
by the notion of a deformation vector which is the solution of the system db/da = — Mf, 
|l30|| . However, for sake of simplicity, we shall consider only the cases when detM^ 7^ 
almost everywhere in J. The Jordan triple system with this property are called Jordan 
triple system with invertible elements. 

The following Lemma demonstrates that some properties of are the same as in a 
ring. 

Lemma 14. Let the operator Ma be invertible then 



{abc} = {cba}, 

{ab{cde}} — {cd{abe}} = {{cba}de} — {c{bad}e}. 



(4.1) 
(4.2) 



Ma-. = 



a 




1 



= a. 



a 




(4.5) 



and the following rules of cancellation of a and a are valid 



{aa ^5} = 6, {a ^{acb}a ^} = {cba ^}, 
{a{a-^ca-^}b} = {ca-H}, {ba{a-^ ca-^]} = {bca-^}. 



(4.6) 



Proof. Let b e J, c = Ma{b), d = Mq-i(c). One obtains using ( ^ ) 

Ma{d) = {a{a''^ca^^}a} = 2{aa"^{aa"^c}} - {a{a''^aa^^}c} 



= 2{aa {aa {aba}}} — {{aa a}a c} = {aa {aba}} 



= {aba} = Maib), 



and therefore b = d. Arbitrariness of b implies M^-iMq = /. 
Further on, {a^^)~'^ = M-\(a-i) = Ma(a-i) = a. 
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For arbitrary b 

{aa-^b} = {aa-'^{aM-\b)a}} = {{aa-^a}M-^{b)a} = {aM-^{b)a} = b. 
Taking this into account when differentiating relation a = {aa^^a} one proves the last 



formula in (4.5). 



In virtue of (4^) and relation L^^-i = I which is already proved one obtains 



and the last formula in ( |4.6| ) is proved. The next to the last is equivalent (use operator 
notation). Finally, ( [4.2| ) implies 

{a^^ {acb^a^^^} = 2{a^^bc} — _i(6)a(c) = {a^'^bc}. ■ 

Now we can prove the formula ("harmonic mean") 

{a-^ + b-^)-^ = {a{a + b)~^b]. (4.7) 

Let us denote a + 6 = c""*^, then 

M(„-i+f,-i)({ac6}) = {a^^ {acb}a'^} + {b^^ {acb}b'^} 
+ 2{a-^{ac[c-^ - a)}b-^} 
= {cba~^} - {r^ac} + 26-1. 

Symmetrization on a and b gives (|4.7D . 

Analogue of the Killing form in the Jordan triple system is the scalar product (a, b) = 
tr Lab- Relation ( [4. 41 ) implies the invariance property 

{{abc},d) = {a,{bcd}) (4.8) 

of this product. Further on we assume that it is also symmetric and nondegenerate. Notice 
that if the element a is invertible then the equalities Lab = MaMa-ib, Lba = Afa-ife-^a 
imply (o, 6) = {b,a). Hence, in virtue of continuity, in the Jordan triple systems with 
invertible elements the symmetry property always holds. 

Examples below together with the reductions a = =ba^ of the Examples 1 and 2 (with 
M = N) exhaust all simple Jordan triple systems aside from two exceptional ones. 



Example 1. The linear space J of N x N matrices becomes the Jordan triple system 
with respect to the triple product defined by means of the standard matrix multiplication 
as follows 

{abc} = -{abc + cba). 

The operator Ma is invertible iff det a ^ 0, that is almost everywhere. The element 
coincides with inverse matrix. The subspaces of symmetric and skewsymmetric matrices 
are Jordan triple systems as well. However, in the Jordan triple system of the skewsym- 
metric matrices of odd order the operator Ma is degenerate for all a. The scalar product 
is (a, b) = tr ab. 
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Example 2. Previous example admits generalisation for N x M matrices: 
{abc} = -{ab^c + cb^a), 

where denotes transposition. In particular, if M = 1 then J turns into the A^-dimensional 
vector space with multiplication 

{abc} = ^{{a,b)c+ {c,b)a) 

where (, ) denotes the standard scalar product. However, the operator Ma is not invertible 
for M ^N. 

Example 3. More interesting triple product in the A^-dimensional vector space is given 
by 

{abc} = {a, b)c + (c, b)a — (a, c)b. 

The scalar product in J coincides with the standard one. Operator Ma, its inverse and 
vector are defined by formulae 

Ma{b) = 2{a, b)a - (a, a)6, M^^ = (a, a)~'^Ma, a'^ = {a, ay^a. 
4.2 Jordan analogues of Heisenberg equations 

The lattice ( |2.2D , (A) of discrete relativistic Toda type admits literal generalization 

^{Tm - + HTn - l)y^^ + KTmTn " 1)^"^ =0, \ + ^i + U = {) (4.9) 

for arbitrary Jordan triple system with invertible elements. Indeed, using relation ( [4.7D 
one can prove that the mapping given by the formulae 

X = vy-^ + \{x + vY^, Y = -fix'^ - X{x + y)-^ 

coincide with its inverse and therefore the duality transformations ( ^.4| ) map equation 
([4.9| ) into itself. The nonrelativistic analogue of this equation corresponding to the set of 
parameters /x = 1, v = — 1, A = can be obtained along the same scheme as in the Section 



2.3 



Equation (^^) is the Euler equation for the Lagrangian 
C = ^{fi log det M^ + u log det My + X log det M^). 

m,n 

In order to prove this, note that if f{u) = ^ log det then, in virtue of Lemma 14, 

(^,1^) = -^f{u + £v)\e=o = ti:{M~^Muv) = trL„-i^ = {u~^,v) 
and therefore df /du = u~^. 
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Anafogously, the lattice (fo) of the relativistic Toda type admits generalization 
ij = p-p-i, p = Mp{M-\pi) - M-\p^i) - 1 + y-i) 
corresponding to the Lagrangian 

C= j dt log det Mp - {p, y-^) - ^ log det My). 

Indeed, if f{u) = {u^^,a) then 

^fiu + ev)\e=o = -{M-\v),a) = -{M-\a),v) 

and therefore df /du = —M^^{a). 

The duality transformations (3.4) are defined by the mapping T of the form 



Y=p-^-y-\ P = y-'-M-\p) 



which is involutive. Applying the scheme of the Section 3.3 one obtains the Toda type 
lattice 

q = -M^{{q, - g)-^ - (g - q_,y^) 



which is equivalent to the Jordan Volterra lattice |2£] 
p = MpiP-Pi), P = Mp{p_i-p). 
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